Nuclear symmetry energy probed by neutron skin thickness of nuclei 
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We describe a relation between the symmetry energy coefficients c sym (p) of nuclear matter and 
a S ym(j4) of finite nuclei that accommodates other correlations of nuclear properties with the low- 
density behavior of c sym (p). Here we take advantage of this relation to explore the prospects for 
constraining c aym (p) of systematic measurements of neutron skin sizes across the mass table, using 
as example present data from antiprotonic atoms. The found constraints from neutron skins are in 
harmony with the recent determinations from reactions and giant resonances. 

PACS numbers: 21.10.Gv; 21.65. Ef; 21.60.-n; 26.60.-n 



A wealth of measured data on densities, masses and 
collective excitations of nuclei has allowed to resolve ba- 
sic features of the equation of state (EOS) of nuclear 
matter, like the density po w 0.16 fm~ 3 , energy per par- 
ticle a v « —16 MeV, and incompressibility K v w 230 
MeV |l| at saturation. However, the symmetry prop- 
erties of the EOS due to differing neutron and proton 
numbers remain more elusive to date. The quintessen- 
tial paradigm is the density dependence of the symmetry 
energy 0, 0, 0, 0, 0, 0, 0, Qiij] . The accurate charac- 
terization of this property entails profound consequences 
in studying the neutron distribution in stable and exotic 
nuclei and neutron-rich matter [3, 0, 0]. It impacts on 
heavy ion reactions [B|, 0, 0, 0, 0], nuclear astrophysics 

0, El ! and on diverse areas such as tests of the Stan- 
dard Model via atomic parity violation [ll[. 

The general expression e(p, 5) = e(p, 0) + c sym (p)<5 2 + 
0(S 4 ) for the energy per particle of nuclear matter of 
density p — p n + p p and asymmetry 5 — (p n — p p )/p de- 
fines the symmetry energy coefficient c sy m(p) of a nuclear 
EOS. It is customary and insightful to characterize the 
behavior of an EOS around the saturation density po m 



terms of a few bulk parameters, like e(p, 0) 
and c. 
5. 6. 



~ J-Le- 
The value of J 



\K sym e 2 where e = (p -p)/(3p ) 



(p ) is acknowledged 



to be about 32 MeV. The values of L — 3pdc sym (p) /dp\ 



constraint c sym (0.1fm 6 ) = 23.3-24.9 MeV [1||, imply- 
ing 7 ~ 0.5-0.65. Note that the Thomas-Fermi model fit- 
ted very precisely to binding energies of 1654 nuclei [13] 
predicts an EOS that yields 7 = 0.51. With the caveat 
that the connection of experiments to the EOS often is 
not at all trivial [1, 0, 0, 0, El , it is important to seek 
further clues from the above and other isospin-sensitive 
signals, such as the neutron skin thickness S = R n — R p 
of nuclei (difference of neutron and proton rms radii). 
Because S of heavy nuclei correlates linearly with the 
slope L of Csymin mean field theories of nuclear structure 
0, 0, 0, 0, EH El , these studies have far-reaching 
implications for nuclear theory. 

In this work we show that c sym (p) of the EOS equals at 
p w 0.1 fm~ 3 the value of the symmetry energy coefficient 
asym(^4) of heavy finite nuclei, universally in mean field 
theories. The observed correlations of S 0,0,0, 0, 0, 
and of the excitation energy of the GDR [lj] with the 
density dependence of c sym can be deduced naturally 
from this relation. We resort to the nuclear droplet model 
(DM) [I2J to work out the analytical formulas. The re- 
sult derived for S is applied to investigate limits to the 
£ 2 slope and curvature of c sym from neutron skins measured 
for 26 stable nuclei, from 40 Ca to 238 U, in antiprotonic 
atoms 20]. A main point is ascertaining how far uni- 
formly measured neutron skins over the periodic table 



and K synl = 9p 2 d 2 c sym (p) /dp 2 \ po govern the'density de- ^ . hel P constrain the density dependence of c sym . We 
pendence of c sym around po- They are less certain and 
the predictions vary largely among nuclear theories, see 
e.g. Ref. [3] for a review. 

In experiment, recent research in intermediate-energy 
heavy ion collisions (HIC) is consistent with a depen- 



dence c sym (p) - c sym ( Po )(p/po) 7 at p < po @,0,ii. 
Isospin diffusion predicts 7 = 0.7-1.05 (L = 88 ± 25 
MeV) @, 0], isoscaling favors 7 = 0.69 (L - 65 MeV) 
Q, and a value closer to 0.55 (L ~ 55 MeV) is in- 
ferred from nucleoli emission ratios 0]. Nuclear reso- 
nances are another hopeful tool to calibrate c sym (p) be- 
low p [3, 14 , HI, III . Indeed, the giant dipole resonance 
(GDR) of 208 Pb analyzed with Skyrme forces suggests a 



provide first evidence that the constraints from skins are 
in consonance with the recent observations from reac- 
tions and giant resonances, though the probed densities 
and energies are not necessarily the same. 

The symmetry energy coefficient a sym (^4) of finite nu- 
clei is smaller than the bulk value J. Given a nuclear 
force, the DM allows one to extract a sym (A) as [12l. |21| 



J 



1 + XA 



with xa 



(1) 



The so-called surface stiffness Q measures the resistance 
of the nucleus against separation of neutrons from pro- 
tons to form a neutron skin. One can obtain Q of nu- 
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FIG. 1: (Color online) Correlation of the quantal selfconsis- 
tent S value in 208 Pb with the slope of the symmetry energy 
L (a), the ratio L/J (b), and with J — a sym (A) (c), for various 
nuclear models (DD and PC stand for density dependent and 
point coupling models). From left to right, the correlation 
factors are r = 0.961, 0.945 and 0.970. 



clear forces by asymmetric semi-infinite nuclear matter 



(ASINM) calculations [12], [21 



|. The contribution of 
a sym (A) to the nucleus energy is a sym (A) (I + xa!c) 2 A, 



^sym V 
o2 ' 



where I = (N - Z)/A and I c = e 2 Z/(20JR) is due to 
Coulomb. One has R = r^A 1 / 3 . A small correction to 
asym(^l) from surface compression [l2| is neglected here. 
Let us mention that ((T|) may be derived also from the 
notion of surface symmetry energy [1, . 

The neutron skin thickness of nuclei is obtained as 



2 Z/(70J)] + S s 



(2) 



in the DM [12], |23[. The quantity t gives the distance 
between the neutron and proton mean surface locations: 



t 



3r J/Q 



2 1 + x A 
2ro 
3J 



[J- 



(I-Ic) 
sym (A)} A 1 / 3 (I - Ic), 



(3) 



where in the second line we have introduced the surface 
symmetry term a SB (A) = [J — a sym (A)]A 1 / 3 using Eq. fl}. 
The second term in Eq. ([2|) is due to Coulomb repulsion, 
and S^w = ^3/5 [5(6 2 — b 2 ) / (2i?)] is a correction caused 
by an eventual difference in the surface widths b n and b p 
of the neutron and proton density profiles. 

We first illustrate the aforesaid correlation of S of 
heavy nuclei with L in Fig. 1(a). It depicts the quantal 
self-consistent value of S in 208 Pb against L for multiple 
Skyrme, Gogny, and covariant models of different nature 
@, H i, H, 0, 0, 0, EH, 0). In Fig. 1(b) we show that 
a similar correlation exists with the ratio L/J, which is 
proportional to 7 if a scaling (p/po) 7 holds for c sym (p). 
And in Fig. 1(c) we show that the close dependence of S 
on J — a sym (A) predicted by the DM is borne out in the 
quantal S value, using forces where we have computed 
Q in ASINM. It reassures one that the DM expression 



TABLE I: Value of J, &sym (^4) &nd. density p that fulfils 
Csym (p) = a sym (A) for A = 208, 116 and 40, 



m various nu- 
3 . Here 



clear models. J and a sym are in MeV and p is in fm 
Csym(p) was computed exactly as ^ d 2 e(p, S)/dS 2 \s=o from the 
EOS of the models. 
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0.105 


24.6 


0.099 


21.3 
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FSUGold 


32.6 


25.4 


0.099 


24.2 


0.092 


21.9 


0.078 


TF [17] 


32.6 


24.2 


0.094 


22.9 


0.086 


20.3 


0.071 


SLy4 


32.0 


25.3 


0.100 


24.2 


0.093 


22.0 


0.079 


SkX 


31.1 


25.7 


0.103 


24.8 


0.096 


22.8 


0.084 


SkM* 


30.0 


23.2 


0.101 


22.0 


0.094 


19.9 


0.079 


SIII 


28.2 


24.1 


0.093 


23.4 


0.088 


21.8 


0.078 


sen 


26.8 


21.6 


0.104 


20.7 


0.098 


18.9 


0.084 



incorporates the proper elements for the study. Many 
of the given nuclear interactions are accurately fitted to 
experimental binding energies, single-particle data, and 
charge radii of a variety of nuclei. However, their isospin 
structure is not sufficiently firmed up as seen e.g. in the 
differing predictions for S'( 208 Pb). There is thus a need 
to deepen our knowledge of isospin-sensitive observables 
like S and of their constraints on c S y m (p). 

We bring into notice a genuine relation between the 
symmetry energy coefficients of the EOS and of nuclei: 
c S ym(p) equals a syin (A) of a heavy nucleus like 208 Pb at 
a density p ~ 0.1 fm~ 3 . Indeed, the relation holds sim- 
ilarly down to medium mass numbers, at lower p values 
and a little more spread. Table I exemplifies this fact 
with several nuclear models, where we show the density 



fulfilling c sym {p) = a sym (A) for A = 208, 116, and 40. 
We find that this density can be parametrized as 

PA = Po-po/(l + cA 1 / 3 ) (4) 

with c fixed by P208 = 0.1 fm -3 (which gives p n6 ~ 0.093 
fm~ 3 and p 40 ~ 0.08 fm~ 3 for the models of Table I). 



> 40 ~ 

The relation 



»(P) 



i(A) n can be very help- 



ful to elucidate other correlations of isospin observables 
with c sym (/?) and to gain deeper insights into them. For 
example, it allows one to replace a sym (A) in Eq. ([3]) for 
a heavy nucleus by c sym (p) ~ J — Le 
computed at p ~ 0.1 fm~ 3 [25]: 



\K sym e 2 with e 



t = 



2r / _ ^sym 

3J V 6 2L 



lA^d-Ir 



(5) 



The imprint of the density content of the symmetry en- 
ergy on the neutron skin appears now explicitly. The 
leading proportionality of (O with L explains the ob- 
served linearity of S of a heavy nucleus with L in nuclear 
models @, 0, 0] ■ The correction with K sym does not alter 
the situation as e ~ 1/9 is small. One can use Eq. §5§ in 



3 



other mass regions by calculating e from pa of Eq. ([4}. 
We have checked numerically in multiple forces that the 
results closely agree with Eq. © for the 40 < A < 238 
stable nuclei given in Fig. 2. 

With the help of Eq. (|5|) for t (using pa to compute e), 
we next analyze constraints on the density dependence 
of the symmetry energy by optimization of @ to exper- 
imental S data. We employ c sym (p) = 31.6(/o/po) 7 MeV 
S IH El and take as experimental baseline the neutron 
skins measured in 26 antiprotonic atoms (see Fig. 2). 
These data constitute the largest set of uniformly mea- 
sured neutron skins over the mass table till date. With 
allowance for the error bars, they are fitted linearly by 
S = (0.9 + 0.15)7+ (-0.03 + 0.02) fm 0. This systemat- 
ics renders comparisons of skin data with DM formulas, 
which by construction average the microscopic shell ef- 
fect, more meaningful (2(|. We first set b n — b p (i.e., 
S S w = 0) as done in the DM 12, 2^, 26[ and in the anal- 
ysis of data in Ref. [l^. Following the above, we find 
L = 75 ± 25 MeV (7 = 0.79 ± 0.25). The range AL = 25 
MeV stems from the window of the linear averages of 
experiment. The L value and its uncertainty obtained 
from neutron skins with S sw = is thus quite compat- 
ible with the quoted constraints from isospin diffusion 
and isoscaling observables in HIC @, 0, S] ■ On the other 
hand, the symmetry term of the incompressibility of the 
nuclear EOS around equilibrium (K = K V +K T 8 2 ) can be 
estimated using information of the symmetry energy as 
K T « K sym -6L The constraint K T = -500±50 



MeV is found from isospin diffusion @, 0] , whe 



reas our 



study of neutron skins leads to K T = — 500^ioo MeV. A 
value K T = —550 ± 100 MeV seems to be favored by the 
giant monopole resonance (GMR) measured in Sn iso- 
topes as is described in [l3| ■ Even if the present analyses 
may not be called definitive, significant consistency arises 
among the values extracted for L and K T from seemingly 
unrelated sets of data from reactions, ground-states of 
nuclei, and collective excitations. 

To assess the influence of the correction S sw in ( ffl) w e 
compute the surface widths b n and b p in ASINM [22j. 
This yields the 6„( p ) values of a finite nucleus if we re- 
late the asymmetry 60 in the bulk of ASINM to I by 
S a (l + xa) = I + xaIc 21, 22, 23|]. In doing so, we find 



that Eq. ([2|) reproduces trustingly S (and its change with 
/) of self-consistent Thomas-Fermi calculations of finite 
nuclei made with the same nuclear force. Also, S sw is 



very well fitted by S s . 



(J SW I. All slopes a sw of the 



forces of Fig. 1(c) lie between cr™ n = 0.15 fm (SGII) and 
a sw X = 0-31 f m (NL3). We then reanalyze the exper- 
imental neutron skins including S™ n and S™ K in Eq. 
([2]) to simulate the two conceivable extremes of S sw ac- 
cording to mean field models. The results are shown in 
Fig. 3. Our above estimates of L and K T could be shifted 
by up to —25 and +125 MeV, respectively, by nonzero 
5 SW . This is on the soft side of the HIC [1,0, @| and 
GMR [l3j analyses of the symmetry energy, but closer 
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FIG. 2: (Color online) Comparison of the fit described in 
the text of Eq. ([2]) with the experimental neutron skins 
from antiprotonic measurements and their linear average S = 
(0.9 ± 0.15)1 + (-0.03 ± 0.02) fm 0. Results of the modern 
Skyrme SLy4 and relativistic FSUGold forces are also shown. 
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FIG. 3: (Color online) Constraints on L and K T from neutron 
skins and their dependence on the S BV , correction of Eq. ([2| . 
The crosses express the L and K T ranges compatible with the 
uncertainties in the skin data. The shaded regions depict the 
constraints on L and K T from isospin diffusion [fj [?J and on 
K T as determined in [l3|] from the GMR of Sn isotopes. 



to the alluded predictions from nucleon emission ratios 
0, the GDR 14 1, and nuclear binding systematics [l7j |. 
One should mention that the properties of c sy m(p) de- 
rived from terrestrial nuclei have intimate connections to 
astrophysics [1, As an example, we can estimate 

the transition density pt between the crust and the core of 
a neutron star @,[n3| as pt/po ~ 2/3 + (2/3)^ K xvm /2K V , 
following the model of Sect. 5.1 of Ref. [lfj. The con- 
straints from neutron skins hereby yield p t ~ 0.095 + 0.01 
fm~ 3 . This value would not support the direct URCA 
process of cooling of a neutron star that requires a higher 
Pt [EE]- The result is in accord with p t ~ 0.096 fm -3 
of the microscopic EOS of Friedman and Pandharipande 
[27| |. as well as with p t ~ 0.09 fm -3 predicted by a recent 
analysis of pygmy dipole resonances in nuclei [151 ] . 

We would like to close with a brief comment regard- 
ing the GDR. As mentioned, Ref. 1J] very interestingly 
constrains c sym (0.1) from the GDR of 208 Pb. The anal- 



4 



ysis notes that the mean excitation energy of the GDR 
depends on g(A) = J /{I + |a ss (A)A -1 / 3 / J} [J[l| and 
shows numerically that the values of <?(208) and c sym (0.1) 
are correlated in Skyrme forces. Inserting a ss (j4) given 
below Eq. ©, one has g{A) = J/{1 + ^[J -a sym (A)]/J}. 
Immediately, the equivalence a sym (208) ~ c sym (0.1) ex- 
plains why g(208) has a dependence on c sym (0.1), gives 
it anal ytic ally, and validates it for any type of mean field 
model |28| . One could extend it to other A values through 
Eq. ((4]). In conclusion, the discussed relation of c sym (p) 
with a sym (^4) can be much valuable to link different prob- 
lems depending upon a sym (A) of nuclei to the symmetry 
properties of the EOS. 

Summarizing, we have described a generic relation be- 
tween the symmetry energy in finite nuclei and in nuclear 
matter at subsaturation. It plausibly encompasses other 
prime correlations of nuclear observables with the den- 
sity content of the symmetry energy. We take advantage 
of this relation to explore constraints on c sym (£) from 
neutron skins measured in antiprotonic atoms [20J. We 
discuss the L and K T values that skins favor vis-a-vis 
most recent observations from reactions and giant reso- 
nances. The difficult experimental extraction of neutron 
skins limits their potential to constrain c sym (p). Interest- 
ingly, we learn that in spite of present error bars in the 
data of [20j , the size of the final uncertainties in L or K T 
is comparable to the other analyses. This highlights the 
value of having skin data consistently measured across 
the mass table, and calls for pursuing extended mea- 
surements of neutron radii and skins with "conventional" 
hadronic probes. Combined with a precision extraction of 
R n of 208 Pb through electroweak probes [29(, they would 
contribute to cast uniquely tight constraints on c sym (p). 
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